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COMPUTER  SIMULATION  OF  INTENSE  ELECTRON  BEAM 
GENERATION  IN  A  HOLLOW  CATHODE  DIODE 


Introduction 

Several  years  ago  at  the  Naval  Research  Laboratory* Cooper stein,  et  al., 
observed  some  interesting  anode  damage  patterns  during  runs  of  the  Gamble  I 
low-imped2mce  diode. ^  The  experimental  arrangement  is  depicted  in  Fig.  1. 

A  simple  carbon  hollow  cathode  opposed  a  flat  plate  porous  graphite  anode. 
(Note  that  only  the  small  D  carbon  witness  plate  location  is  of  interest 
here. )  Large  external  magnetic  field  coils  maintained  the  entire  configura¬ 
tion  in  a  30  kG  solenoidal  magnetic  field.  A  number  of  discharges  were  nin 
in  the  600  kV  to  1  MV  range.  After  a  typical  experiment,  exeunination  re¬ 
vealed  that  a  perfect  circular  groove  had  been  cut  into  the  anode  surface. 

The  dimensions  of  the  groove  were  such  that  the  cathode  could  neatly  fit  into 
it.  Furthermore,  upon  closer  examination  there  was  an  even  finer  pattern  of 
damage  inside  the  groove.  A  radial  cross-section  had  the  shape  of  a  straights 
sided  "W".  The  groove  was  decidedly  deeper  at  the  edges  them  along  the 
center. 

The  close  correspondence  between  the  cathode  radius  and  the  damage 

radius  was  attributed  to  the  effect  of  the  strong  axial  magnetic  field.  This 

effect  alone,  however,  could  not  explain  the  sheer  magnitude  of  the  damage. 

2 

Some  100  kA/cm  of  electron  current  density  would  have  been  required  to 

pulverize  the  anode  material,  standard  one-dimensional  analysis  cannot 

explain  the  emission  of  such  high  current  densities  from  the  cathode. 
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Nevertheless,  their  existence  has  been  established.  A  possible  explanation 

is  that  electron  emission  is  enhanced  through  a  two-dimensional  intensifica- 
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tion  of  the  electric  field  near  the  cathode  face  edges.  '  This  enhanced 
edge  emission  would  also  explain  the  ”W"  shape  of  the  anode  groove  cross- 
section.  The  computer  simulations  described  below  were  conducted  to  test 
this  theory. 

The  Simulation  Code 

The  computer  code  utilized  in  these  studies  is  a  modified  version  of  the 

DI0DE2D  code  developed  at  the  Naval  Research  Laboratory  during  the  mid-1970's 

by  one  of  the  authors  (R.E.L.).  That  original  code  performed  numerous 

accurate  and  efficient  one-  and  two-dimensional  electrostatic  diode  simula- 
4  5 

tions.  '  These  results  encouraged  the  expansion  of  the  program  to  allow  the 
full  treatment  of  radial  and  axial  self-magnetic  fields  as  well  as  the  inclu¬ 
sion  of  an  imposed  axial  magnetic  field.  As  in  the  original  code,  great 
emphasis  was  placed  upon  optimizing  the  numerics  of  the  new  code  for  use  on 
NRL's  vector  machine,  the  Texas  Instrument  Advanced  Scientific  Computer.  The 
coding  work  was  completed  by  another  of  the  authors  (R.J.B.)  in  August  of 
1978.  Since  then,  work  to  improve  its  efficiency  has  continued  simultaneous 
to  its  application  to  a  variety  of  problems  of  interest  including  the  one 
presented  herein.  The  capabilities  of  the  code  2md  its  limitations  are  out¬ 
lined  below. 

The  code  is  2h  dimensional,  particle-in-cell  (PIC).  Inhomogeneities 
are  allowed  in  the  radial  (r)  ^md  eocial  (z)  spatial  dimensions.  Con^lete 
azimuthal  symmetry  is  assumed.  On  the  other  hand,  all  three  momentum  com¬ 
ponents  (r,  6,  z)  are  retained.  This  retention  of  the  theta  momentum  is  a 
prerequisite  for  the  treatment  of  radial  and  axial  magnetic  fields  in  the 
simulation.  The "particles"  in  this  model  are  axially-centered  rings  of 
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charge  with  all  other  degrees  o£  freedom.  Xn  reality  they  are  macroparticles 
carrying  many  times  elementary  charge  but  retaining  the  physical  charge"- 
to-mass  ratios  of  the  protons  and  electrons  which  they  represent.  Area 
weighting  (i.e.,  linear  interpolation)  is  used  to  couple  these  charges  with 
the  electric  and  magnetic  fields  calculated  over  a  fixed  set  of  grid-points 
in  the  region  of  interest.  The  fields  thus  interpolated  to  the  particle 
positions  act  on  these  charge-current  rings  by  way  of  the  relativistic 
Lorentz  force  equation. 

The  fields  are  treated  electrostatically  emd  magnetostatically.  In 
this  sense  the  code  does  not  perform  a  true  "simulation"  since  Maxwell's 
equations  are  not  observed.  Rather  the  treatment  is  "quasistatic".  Equilib¬ 
rium  solutions  to  various  diode  geometries  are  sought.  The  "timesteps" 
which  appear  in  the  code  are  actually  snapshots  of  the  system  tdiile  it  seelcs 
to  relax  toward  its  steady-state  operation.  It  is  this  steady  state  tdiich 
is  of  primary  interest.  In  order  to  determine  the  electric  field  within  the 
diode  region  particle  charge  densities  are  distributed  over  a  fixed  grid  and 
the  discrete  Poisson's  equation  is  solved.^  The  same  technique  is  used  to 
determine  the  axial  and  radial  components  of  the  self-magnetic  field  from  the 
azimuthal  current  densities.  In  both  cases  the  code  permits  irregular 
conducting  boundaries  inside  the  computational  region.  The  treatment  of  such 
internal  boundaries  entails  the  use  of  a  "capacitance  matrix".^  The  internal 
surfaces  thus  created  are  held  at  predetermined  electric  potential  values 
while  magnetic  flux  from  all  but  the  imposed  fields  is  excluded  from  them. 
Direct  radial  integration  of  the  axial  current  densities  over  the  smsh  yields 
the  azimuthal  magnetic  field  via  Aspere's  Law.  The  outer  radial  boundary  of 
the  diode  region  may  be  either  conducting  or  free-space.  Hie  electrostatic 
potential,  is  set  constant  along  all  conducting  boundaries.  Along  radial 


free-space  boundaries,  is  graded  logarithmically  emd  along  axial  ones, 
linearly.  For  A.,  when  a  free-space  outer  radial  boundary  is  chosen,  the 
external  source-free  cylindrical  vector  potential  eigenvectors  are  matched 
up  with  their  internal,  source-consistent  counterparts  there.  The  algorithm 
also  permits  the  imposition  of  a  temporally  and  spatially  uniform  axial 
magnetic  field.  This  feature  was  essential  to  the  completion  of  this 
particular  study. 

At  the  start  of  a  typical  computer  run,  the  computational  diode  region 
is  a  complete  vacuum,  devoid  of  particles.  The  electric  potential  is  pre-set 
along  the  entire  boundary  as  well  as  along  all  internal  conducting  surfaces. 

A  fixed  value  is  similarly  set  for  the  imposed  axial  magnetic  field  (B^) . 

The  emission  of  the  ions  is  permitted  anywhere  along  the  conducting  high 
voltage  boundary  (anode)  surface.  Electrons  are  emitted  along  the  front  face 
of  the  cathode  surface.  The  value  of  the  perpendicular  electric  field  at  a 
given  emission  point  determines  the  total  clfiarge  (i.e.,  number  of  particles) 
that  will  be  emitted  there.  At  the  beginning  of  a  timestep,  the  electric 
field  at  a  surface  specifies  the  net  charge  density  on  the  surface  via  Gauss* 
Law.  The  surface  integral  of  this  density  over  a  cell  width  around  a  given 
grid  point  yields  the  net  charge  which  is  emitted  there  for  this  timestep. 
Prior  to  the  actual  particle-pushing  the  electric  field  is  recalculated 
talcing  into  account  the  newly  emitted  charge. 

The  emitted  particles  are  then  pushed  according  to  the  relativistic 
Lorentz  force  law  using  the  area  weighted  electric  and  magnetic  field  values 
interpolated  at  the  particle  position  from  the  four  nearest  grid  points. 
After  the  pushing  in  each  timestep  the  charge  and  current  density  associated 
with  each  particle  is  distributed  over  the  four  nearest  gridpoints  using  the 
same  linear  interpolation  scheme  in  reverse.  This  yields  a  complete  array  of 


values  for  the  charge  density,  p,  and  the  current  density,  J,  over  the  con^u- 
tational  mesh.  Poisson  solving  these  arrays  yields  ^  emd  A.  from  which  E-, 

V  R 

E  ,  B  and  B  are  calculated.  The  azimuthal  coiqponent  of  B  is  obtained 

Z  R  Z 

through  direct  integration  of  J  over  the  grid.  Quantities  of  interest  are 

z 

then  extracted  and  output  via  diagnostic  subroutines.  The  code  then  cycles 
to  the  next  tlmestep  for  particle  emission. 

Finally,  it  would  be  noted  that  the  numerics  of  the  particle  pushing  as 
well  as  the  potential  solving  has  been  completely  "vectorized".  Thus  the 
momentum,  position,  and  field  components  associated  with  the  entire  ensemble 
of  particles  are  treated  as  macro-vector  quantities.  Arithmetic  operations 
performed  with  them  are  accomplished  in  a  conipletely  vector-array  format. 

This  property  of  the  code  permits  efficient  running  times  on  the  most 
advanced  scientific  con^uters.  (Of  course,  the  interpolation  of  p,  J,  B,  and 
B  values  between  particle  positions  euid  grid  points  requires  random  accessing 
of  array  points  and  this  process  cannot  be  vectorized.) 

The  Computer  Experiment 

Configuration  Uhder  Study 

The  diode  configuration  modeled  in  this  simulation  was  chosen  specif¬ 
ically  for  its  significant  distortion  of  the  internal  electric  field.  Any 
emission  enhancement  resulting  from  this  distortion  should  be  relatively  easy 
to  see.  Figure  2  provides  a  scale  drawing  of  the  actual  geometry  of  the 
diode  modeled.  The  hollow-cathode  projects  into  the  computational  region 
from  the  left-hand  boundary.  A  cylindrical  anode  cup  con^letely  surrounds  it 
radially  and  to  the  right.  One-quarter  centimeter  separates  the  cathode  face 
from  the  ^ulOde  disc  plate  while  radial  gap  between  these  two  coaxial  elec¬ 
trodes  is  one  and  one-quarter  centimeters  wide.  The  cathode  itself  measures 
1.27  centimeters  in  outer  radius  with  a  wall  thickness  of  0.14  cm.  It 
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projects  0.08  cm  Into  the  2-0  computational  region.  A  conducting  plug  fills 
the  cathode  interior  to  within  0.08  cm  of  the  edge  face  and  thus  coincides 
with  the  computational  region's  inner  left  boundary.  The  upper  half  of  the 
region's  left  boundary  is  left  "open"  with  the  electrostatic  potential 
graded  logarithmically  with  radius. 

The  vacuum  diode  eguipotential  svir faces  are  depicted  in  Fig.  3.  Of 
particular  interest  in  this  plot  is  the  numner  in  which  these  surfaces  are 
compressed  near  the  edges  of  the  cathode  face.  This  bunching  gives  rise  to 
the  electric  field  enhancement  at  these  points  and  to  the  resulting  enhanced 
electron  emission  there. 

Operational  Parameters 

For  all  of  the  runs  discussed  herein,  the  diode  voltage  is  initialised 

at  120  kV  and  climbs  linearly  over  50  timesteps  to  a  steady  state  plateau 
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value  of  600  kV.  The  timestep.  At,  itself  is  fixed  at  1.6  x  lo  second. 
The  two-dimensional  numerical  mesh  over  which  the  particles  move  measures 
128  cells  in  the  radial  dimension  and  64  cells  in  the  axial.  In  addition,  a 
monolayer  of  "guard  cells”  completely  surrounds  the  mesh  to  yield  a  total  of 
130  X  66  cells  overall.  All  boundary  surfaces  run  between  neighboring  rows 
of  gridpoints  rather  than  along  them.  (Gridpoints  are  defined  as  the  cell 
centers.)  The  electrostatic  potentials  on  these  boundaries  are  held  fixed 
relative  to  the  full  diode  potential  at  any  given  timestep  throughout  the 
simulation.  The  cell  widths  (or  gridpoint  spacings)  are  0.0195  and  0.00521 
centimeters  for  AR  amd  AZ,  respectively.  The  front  face  of  the  cathode 
emits  electrons  over  six  radial  cells.  This  emission  profile  provided 
adequate  resolution  for  this  experiment.  Finally,  a  solenoldal  magnetic 
field  is  loosed  and  is  assumed  to  have  been  on  long  enough  to  have 


completely  penetrated  the  anode  and  the  cathode  before  particle  emission  is 
initiated. 

For  the  first  300  timesteps  of  the  simulation,  only  electron  emission 
was  permitted.  By  that  time,  an  electrons-only  steady  state  had  been 
attained.  The  program  diagnostics  were  then  examined  to  determine  \diich 
regions  of  the  anode  surface  were  being  bombarded  by  the  largest  electron 
current  densities.  At  timestep  301,  ion  emission  was  turned  on  over  those 
regions.  (In  the  physical  diode,  these  are  the  areas  where  the  most  heating 

g 

of  the  anode  surface  has  occurred. }  The  computer  simulation  then  continues 
until  the  net  ion  and  electron  currents  equilibrate.  The  numerical  data 
presented  in  the  next  section  are  taken  from  this  final  timestep  with  the 
device  in  its  steady-state  operation. 

Results 

In  order  to  gauge  the  specific  effect  of  a  strong  imposed  solenoidal 

magnetic  field,  B^,  on  diode  operation,  the  same  diode  was  simulated  with 

and  without  B  =  30  kG.  In  both  cases  the  plateau  voltage  is  set  at  600  kV. 
z 

In  Case  1,  B^  is  set  equal  to  zero.  The  major  phenomenum  to  be  observed  in 
this  case  is  the  two-dimensional  electric  field  enhanced  emission  of  electrons. 
The  results  are  graphically  presented  in  Fig.  4.  After  300  timesteps  of 
electrons-only  operation,  the  equilibrium  profiles  of  electron  emission  cur¬ 
rent  density  at  the  cathode  face  and  the  electron  impact  current  density  at 
the  anode  are  plotted  as  shown.  It  is  quite  clear  that  the  compression  of 
electrostatic  equipotential  liner  at  the  cathode  face  edges  (see  Fig.  3)  has 
had  its  effect.  Electron  emission  is  up  by  a  factor  of  about  2.5  at  the 
edges  as  compared  to  the  center  of  the  cathode  face.  A  peak  electron  current 
density  of  over  78  kA/cm^  is  observed  at  the  edge.  The  sharpness  of  this 
profile  degenerates  markedly  before  hitting  the  anode  face.  Virtually  all 


of  the  double-peaked  definition  is  lost.  The  distribution  is  relatively  dif- 

2 

fuse  with  a  maximum  density  of  only  about  24  kA/cm  .  A  net  current  of 

37.6  kA  now  flows  through  the  diode.  Since  the  critical  current  is  about 

82  kA,  only  very  weak  pinching  can  occur. 

At  timestep  301,  ion  emission  begins  on  the  anode  face  in  the  region 

where  the  electrons  have  been  arriving.  The  addition  of  protons  into  the 

system  has  a  substantial  effect  (see  Fig.  5) .  At  equilibrium,  the  total 

current  through  the  diode  has  jumped  to  61.5  kA.  The  electron  emission  pro- 

2 

file  has  about  the  same  shape  but  it  now  exceeds  90  kA/cm  at  the  catl'.ode 

face  edges.  Similarly,  the  anode  electron  impact  plot  is  still  diffuse  but 

2 

it  now  peaks  at  about  32  kA/cm  .  Note  also  that  the  self-pinching  of  the 
electron  beam  is  more  pronounced.  The  mean  radius  of  the  beam  has  decreased 
by  about  one  millimeter  in  its  passage  across  the  two  and  a  half  millimeter 
gap. 

For  the  diode  simulation  Case  2,  a  thirty  kilogauss  solenoidal  magnetic 

field,  B  ,  is  imposed.  The  results  of  the  runs  are  plotted  in  Figs.  6  and  7. 
z 

As  in  the  previous  case,  only  electrons  flow  through  the  system  during  the 

first  300  timesteps.  Figure  6  reveals  once  again  a  double  peaked  emission 

2 

distribution  at  the  cathode  with  a  maximum  value  of  about  96  kA/cm  .  The 

truly  dramatic  change  is  found  in  the  anode  profile.  The  beam's  radial 

diffusion  while  crossing  the  gap  has  been  sharply  restricted  by  the  imposed 

B  .  The  anode  impact  plot  has  about  the  same  shape  as  the  plot  of  the  elec- 
z 

2 

tron  emission.  It  peaks  at  about  83  kA/cm  .  In  this  steady  state  about 

40  kA  is  flowing  through  the  diode.  The  results  become  more  marked  after 

ions  are  permitted  into  the  system  (see  Fig.  7).  Peak  electron  emission 
2  2 

rises  to  108  kA/cm  while  a  maximum  of  92  kA/cm  of  electron  current  density 
is  being  deposited  at  the  anode  face.  (Note  that  the  total  diode  current  is 


For  the  diode  treated  in  this  paper,  V  equals  600  kV  while  d  is  0.25  centi- 

2 

meters.  The  formula  thus  yields  J  ~  17.4  kA/cm  .  This  prediction  comes 

CL 

2 

fairly  close  to  the  20-24  kA/cm  measured  near  the  center  of  the  cathode 

face.  One  would  physically  expect  those  interior  regions  to  most  closely 

resemble  the  one-dimensional  model.  At  the  edges,  on  the  other  hand,  the 

electron  emission  exceeds  the  one-dimensional  J  prediction  by  a  factor  of 

CL 

4  to  5. 
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The  introduction  of  protons  into  the  system  modifies  the  theoretical 

predictions.  This  new  condition  is  termed  one -dimensional  bi-polar  flow. 

The  new  electron  emission  current  density,  is  approximated  by 

BP 


BP 


»  1.86  J 


CL 


2 

Thus,  it  predicts  about  32  kA/cm  for  the  diode  operating  witn  ions  and 

2 

electrons.  This  agrees  quite  well  with  the  32-42  kA/cm  seen  in  the  simula¬ 
tion  for  the  inner  cathode  face  emission.  However,  it  is  off  by  a  factor  of 
2.5  to  3  from  the  edge  emission  figures.  The  two-dimensional  field  enhance¬ 
ment  effects  have  to  drastically  change  the  physics  there  from  the  simple 
1-D  analysis. 

The  second  major  question  concerns  the  effect  of  the  30  kilogauss 
solenoidal  magnetic  field  on  the  electron  current  densities  impacting  the 
anode  face.  Experimental  corroboration  for  the  computational  results  pre¬ 
sented  herein  was  derived  from  the  previously  mentioned  observations  obtained 
from  a  series  of  runs  by  the  Gamble  I  device  at  NRL  in  the  mid  1970 's.  These 
empirical  results  were  presented  by  Cooperstein  to  the  1976  meeting  of  the 
Plasma  Physics  Division  of  the  American  Physical  Society.  Of  those  results, 
one  of  particular  interest  here  is  shown  in  Fig.  8.  Using  a  diode  configura¬ 
tion  closely  approximating  that  modeled  in  this  computer  simulation,  a  double 
peaked  dcunage  pattern  was  observed  on  the  anode  face  at  very  nearly  the 
radius  of  the  hollow  cathode  after  each  discharge.  Such  damage  could  result 
if  the  anode  material  was  being  pulverized  by  the  electron  beam.  The  pulver¬ 
ization  of  carbon  requires  the  deposition  of  2  kcal/g.  Assuming  a  mean 

graphite  density  of  2.5  j/cc  this  translates  to  5  kcal/cc.  The  electron  beam 

2 

in  the  computer  simulation  was  delivering  a  peak  power  of  54  j/cm  -sec  over 

2 

an  area  of  roughly  1.1  cm  .  Thus  the  experimentally  observed  damage  could 
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be  obtained  to  a  depth  of  about  one  millimeter  if  the  simulated  beam  were 


maintained  for  approximately  forty  nanoseconds.  According  to  Fig.  8,  the 
experimental  voltage  and  current  peaks  were  1  MV  and  200  kA,  respectively. 
Comparing  this  to  the  600  kV  and  60  kA  of  the  simulation  would  indicate  that 
the  experimental  diode  could  have  achieved  the  observed  damage  in  twenty 
nanoseconds  or  less. 

In  short,  the  correlation  between  experiment  and  simulation  is  very 
good.  The  empirical  results  can  now  be  better  understood  as  a  consequence 
of  the  combined  effects  of  two-dimensional  electric  field  electron  emission 
enhancement  as  well  as  of  the  electron  path  stiffening  caused  by  the  strong 
axial  magnetic  field. 
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Experimental  arrangement 


ELECTRONS  HITTING  ANODE  CATHODE  ELECTRON  EMISSION 


Bj=0.  ELECTRONS  ONLY 


Fig.  4  -  Steady-ttat*  currant  danalty  for  caaa  I  with  only  alactrona 
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ELECTRONS  HITTING  ANODE  CATHODE  ELECTRON  EMISSION 


Bz  =  0,  IONS  a  ELECTRONS 


R(x  0.2cm) 

VsGOOkV,  |s6l.5kA 

Fig.  S  -  Steady-state  current  density  for  case  1  with  electrons  and  lens 
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ELECTRONS  HITTING  ANODE  CATHODE  ELECTRON  EMISSION 


02  =  30kG,  ELECTRONS  ONLY 


Fig. 


6  -  Stcady-ctct*  current  danaity  for  caaa  2  with  only  alactrona 
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ELECTRONS  HITTING  ANODE  CATHODE  ELECTRON  EMISSION 


Bz  =  30kG,  IONS  ft  ELECTRONS 


R(x  0.2mm) 

V=600kV,  l=59.2kA 


Fig.  7  -  Steady-itate  current  density  for  case  2  with  electrons  and  Ions 
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